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I review Klaus Kinder-Geiger's contributions to the piiysics of relativistic heavy ion 
collisions, in particular, the Parton Cascade Model. Klaus developed this model 
in order to provide a QCD— based description of nucleus-nucleus reactions at high 
^^ ' energies such as they will soon become available at the Brookhaven Relativistic 

fvj ' Heavy Ion Collider. The PCM describes the collision dynamics within the early and 

dense phase of the reaction in terms of the relativistic, probabilistic transport of 
perturbative excitations (partons) of the QCD vacuum. I will present an overview 
K^ , of the current state of the numerical implementations of this model, as well as its 

jj^ ■ predictions for nuclear collisions at RHIC and LHC. 

o 

ps| ' 1 Introduction 

o: 

Q^ . Klaus Kinder-Geiger (KKG), who perished in the crash of Swissair flight 111 

0^ ' near HaUfax on September 2, 1998, was a briUiant theoretical physicist and 

one of my dearest friends. He was one of those few human beings who are 

truly irreplaceable, not because no one else could carry on Klaus' research, 

but because of the unique way he did physics and almost everything else in 

I— { ' life. Klaus combined the abstract mind of the scientist who works on deep and 

r^ , esoteric questions of nature with the wild mind of the artist who is driven to 

create and perform in extraordinary ways. Klaus became famous as a physicist 

for the work he did on the parton cascade model of relativistic nuclear reactions. 

^^ ' But he was equally famous as an extravagantly creative and perceptive person 

in the communities where he lived and made friends. I always thought that 

knowing Klaus was the closest I would ever get to meeting in person one of 

the French existentialists whose novels I had read as a student. 

As many of you know, Klaus very much liked to draw sketches of him- 
self and the world around him. The one entitled "The Wunderbar World of 
KKG" , featured prominently in his home page on the World Wide Web (Figure 
n[). It captures nicely how Klaus viewed himself: always inspired to explore 
unfamiliar territory, uncover new ideas, and gather new experiences. Klaus 
believed that what really counted in life was the unusual, and he dared to live 
by his belief. When he was still a young student at the University of Frankfurt, 
before I knew him, Klaus had tried his hand at painting. His paintings are 
highly expressive and almost haunting; having seen them once one is not likely 
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Figure 1: The "Wunderbar World 
- J of KKG" sketch greeting visitors 

to Klaus' World Wide Web site at 
A:«rS- i^ BNL. 



to ever forget them. One of his paintings, shown in Figure y, is particularly 
remarkable: it anticipates Klaus' greatest contribution to physics. Clearly, the 
subject of the painting must be a parton cascade, as Klaus imagined it in his 
artistic mind in 1986.Q 

The hustler, I am convinced, is depicting Klaus himself. (It is anyone's 
guess whom the two other figures represent. Some similarities with well known 
members of the heavy ion physics community are unmistakable.) Lest anyone 
doubts this interpretation of the painting, let me point out that it is signed 
"CucuruUo '86". Vincent de Cucurullo, short: Vinnie, was Klaus' artistic 
pseudonym by which he signed his paintings. The name of his celebrated 
parton cascade model code, VNI, is a rendering of the artist's name under the 
cloak of a scientific acronym. 



2 The Parton Cascade Model 



The parton cascade model (PCM) was proposed by Klaus and mel^ in 1990-91 







and developed much f 
able series of paperst 



pr during the following years by Klaus in a remark- 
Our aim was to describe the energy deposition, 
thermalization, and chemical equilibration of matter in high energy nuclear 
collisions, and to provide a full space-time picture of the collision up to the 
moment when individual hadrons are formed. The model was, at least orig- 
inally, not conceived as an "event generator" that would predict a full set of 



"^Of course, no one - not even Klaus himself - would have had any notion of this at the time 
the scene was painted. But the artist sometimes envisions novel ideas and concepts long 
before the scientist. What is unusual here is that the artist became the scientist! 




Figure 2: The Pool Players. By Vincent De CucuruUo, 1986 



hadron momentum distributions in the final state. 

In order to enable experimental predictions, Klaus decided to develop 
the parton cascade model code VNI, which contains the implementation of a 
hadronization scheme in the framework of a parton coalenscence model. This 
required a number of compromises and ran somewhat counter to the original 
purpose of the PCM, namely, to explore the range of validity of pcrturbative 
QCD in nuclear reactions. Nevertheless, the predictions of the PCM with an 
added hadronization stage have been, and continue to be, very useful. We 
simply do not know how to do better at the present time. 

The conceptual basis of the PCM is the inside-outside cascade modelH of 
high-energy hadron reactions, which implements the concept that new matter 
produced in hadronic interactions at high energy is formed outside the inter- 

rcting world-tubes of the colliding hadrons. Bjorken's hydrodynamical model 
was devised to describe the evolution of this newly formed matter in the 
central space-time region after thermal equilibration. Beginning in the mid- 
1980's it was realized that the deposition of energy into this region may be, 
at least partially, described in terms of concepts based on perturbatiVje-QCLD 
(minijets) when the energy of colliding heavy nuclei becomes very largepEi'r'^trn 



A computer code (HIJIMG^ incorporating some of these ideas was developed 
by Gyulassy and WangE3lij 

The parton cascade model combined these ideas into one unified scheme 
for the description of the space-time evolution of matter in nuclear reactions. 
Its three main ingredients are: 

1. The initial state is viewed as incoherent ensemble of partons determined 
by the nuclear parton distribution functions qf{x, Q^) and g{x,Q^), where 
the subscript /denotes the quark flavor and g(x,Q^) stands for the gluon 
distribution, x = Pz/ P is the longitudinal momentum fraction of the nu- 
cleon carried by the parton, and Q^ is the parton "scale" or virtuality. 
Before any interaction occurs, Q^ is generally taken as space-like. Our 
knowledge about the space-time structure of the nuclei before the colli- 
sion and our limited information about the intrinsic transverse momenta 
of partons is then used to construct a model for the six-dimensional phase 
space distributions of partons before the interaction: qf (r, p) , g{r^ p) . The 
parton distributions are conveniently initialized at the scale Qq = {p'^)co\i 
of the average momentum scale of the primary parton-parton interac- 
tions. 

2. The time evolution of the parton phase distributions is governed by a rel- 
ativistic Boltzmann equation with "leading-log" improved lowest-order 
collision terms. Only binary interactions are allowed, but the final state 
can have (and generally has) more than two particles. As is well known, 
the higher-order improvement of the cross sections by means of the lead- 
ing logarithmic approximation is equivalent to the scale evolution of the 
parton distributions according to the DGLAP equation. Motivated by 
quantum mechanical considerations, the space-time picture of parton 
propagation before and after interactions is closely related to their off- 
shell propagation: the formation of a parton with virtuality scale Q takes 
a time Tf{Q) « h/Q. 

3. When the parton distributions become sufficiently dilute, they hadronize. 
In VNI the hadronization is described by a clustering algorithm, followed 
by the decay of excited hadrons. The transition is assumed to occur when 
the average virtuality of the partons falls below a critical value Qcrit ~ 1 
GeV, because partons no longer scatter with sufficient energy to allow 
for the collisions to be described by perturbative QCD. 

From a gradient expansion of the evolution equation for the parton Wigner 
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distribution two equations can be derivedllj The first equation 



P'^F^i^'P)-C^i^'P) 



(1) 



describes the free propagation of partons which is intermittently modified by 
interactions given by the binary coUision terms C. The second equation 



P 



_d_ 



Fi{r,p) = Si{r,p) 



(2) 



describes the evolution of the parton distributions with respect to virtuality 
or "off-shellness" p^. This equation is a generalization of the usual mass- 



shell condition F{r,p) ^ 5{p 



to the case where the on-shell particle 



distribution cannot be defined. Si{r,p) describes the splitting of single off-shell 
partons into two partons of smaller virtuality. 

The two equations can be viewed as quantitative representation of Feyn- 
man diagrams of the type shown in Fig. 0. The collision term Ct is represented 
by the binary collision diagram contained in the box at the center of the com- 
plex Feynman diagram, whereas the splitting terms Si are represented by the 
branchings of the initial- and final-state partons .H The Feynman diagram of 
Fig. is finite only if the virtualities of all final-state partons are limited 
from below by some infrared cut-off /i^. In an isolated event ^^ describes 
the hadronization scale, i.e. the virtuality scale below which partons can no 
longer be considered as approximately free, perturbation quanta. In a dense 
medium, where partons rescatter often, p^ is determined by the frequency of 
rescatterings (see Section 3.2). 




Ql 



Figure 3: Graphical representation of the 
QCD transport equations (1,2) defining the 
parton cascade model. The 2 — > 2 scatter- 
ing process at momentum scale Q^ is fol- 
lowed by the virtuality evolution from Q^ 
to /i^. 



In the leading- logarithmic approximation (LLA) , the differential cross sec- 
tion described by the Feynman diagram of Fig. H factorizes into a product of 

''Strictly speaking, Si describes the differential branching probability for an infinitesimal 
change in virtuality; the diagram is an integral representation of Si. 



terms for each of the two incoming and outgoing branch processes and one for 
the binary scattering process. Each branching term, in turn, is represented 
by a product of factors describing the individual branching events and the 
probabiUties for the partons not to branch further in between. In other words, 
the Feynman diagram of the type shown in Fig. y| defines a Markov process 
and, hence, can be described by a probabihstic one-body transport equation. 
This statement is no longer valid, if one tries to go beyond the LLA. However, 
certain effects beyond the LLA can still be described in terms of conditional 
probabilities, such as angular or /c^-ordering of gluons and certain soft-gluon 
interference effects. These effects are quantitatively important and have been 
incorporated into parton cascade codesp 

It is important to note that parton splittings are related by unitarity 
to loop diagrams that describe the running of the strong coupling constant 
cts{Q^)- Both splittings and cts^running are described consistently in the LLA, 
which therefore satisfies the unitarity condition. In plain terms, the combined 
probability for all 2 — > n parton diagrams with n > 3 reduces the probability 
for the occurrence of a 2 — > 2 scattering and so on. By summing all 2 — > n free 
diagrams, but not including the associated loop diagrams, unitarity would be 
violated. This violation leads to the divergence of the sum over n already at 
rather small parton center-of-mass energies, even in the presence of an infrared 
cut-off for the internal propagatorst3 

In the following sections I will discuss two important issues: 

1. The space-time picture governing the initial-state parton distributions. 
This issue is closely connected with the problem of the decoherence of 
the initial parton wavefunctions. 

2. The problem of infrared divergences of the perturbative parton cross 
sections. This issue is closely related to in-medium corrections of these 
cross sections, as well as to coherence properties of the initial-state wave 
functions. 

The approach to local thermal equilibrium has been, extensively studied within 
the framework of the parton cascade picturctd'EJO Without repeating the 
detailed arguments here, let me just state that the PCM approach predicts a 
very short kinetic equilibration time, rth ^ 1 fm/c, which is confirmed by full 
numerical calculationsu 

3 Initial state space— time picture 

The probabilistic interpretation of parton distributions measured in deep- 
inelastic scattering is based on a summation over all final hadronic states. 



A similar interpretation of one-body distributions arising in transport theory 
is grounded on the low-order truncation of the (BBGKY) hierarchy of Green 
functions and on an expansion in powers of h. The validity of this picture ulti- 
mately relies on the separation of time scales in a dynamic process. Although 
these issues are generally well known^^l, their implications for nuclear parton 
cascades have not been fully explored. Recent advances 113 in our understand- 
ing of multiple scattering in QCD have shed some light on the intricacies of 
the formation time concept in non-abelian gauge theories, but it needs to be 
better understood how these results can be consistently incorporated into a 
probabilistic transport theory. 

The original parton cascade model relied, on some basic assumptions about 
initial parton distributions in space-timecJ Denoting the parton light-cone 
momentum by p"*", the parton distributions were assumed to be distributed 
longitudinally according to the uncertainty relation: Ap'^ Ax~ > h. Soft par- 
tons have the widest distributions in the variable x~ , and are assumed to travel 
both ahead of and behind the Lorentz contracted valence quark distributions. 
The argument is that this will not violate causality, because soft partons are 
emitted at a long distance before the collision and, travelling at the speed of 
light, can arrive significantly ahead of the quarks that emitted them. 

The space-time picture of soft partons has been put on a much firmer 
foundation in recent years by the work of McLerr^, Venugopalan and others 
on the random light-cone source model (RLSM)E3 In this model, one views 
the valence quarks constituting the fast-moving nucleus as a thin, Lorentz 
contracted sheet of locally random color sources. The color source is locally 
random, because valence quarks from several nucleons contribute at the same 
point in transverse space. The area density of color sources is given by fi^ = 
Sg^A/irR^, where A is the nucleon number, R the nuclear radius, and g the 
QCD coupling constant. Clearly ^i grows as A^^^, hence can be considered as 
a (potentially) large scale for sufhciently heavy nuclei.^ ^^(/x^) <C 1 can serve 
as the coupling parameter for a new type of perturbative expansion. Formally, 
the model maps into the problem of weakly coupled QCD in the presence of a 
random two-dimensional color source. 

As shown by Kovchegovpl this model can be rigorously derived with stan- 
dard light-cone techniques, which permit an explicit representation of the 
Gaussian ensemble of color sources. This representation can also be used to 
calculate the perturbative emission of soft gluons in collisions between laoxnu- 
clei, described as encounter of two sheet-like clouds of valence quarksEjc3 At 



'^ In practice, fi < 1 GeV even for the heaviest nuclei, even if the "hard" component of the 
gluon distribution is included in the color source. 



leading order this soft gluon radiation is given by: 

dNg _ 4a3 N^^l,^ ^^^^ [ ^^^__F{qa)F{\k~g\a) 



dyd^k^d^h- n^kl N, {TAB{b))J dg ^2(fc_^)2 (3) 

where F{qa) is the color-dipole form factor of_Lhe nucleon and TAsib) denotes 
the nuclear profile function. It can be shown Efl that this classical gluon radi- 
ation matches smoothly onto the perturbative minijet production of gluons at 
higher k^. 

Going beyond the classical approximation by including gluon-loop dia- 
grams leads to a better and more rigorous understanding of the space-time 
distribution of soft gluons in a heavy nucleuses The quantum corrections can 
be formulated in the framework of a space-time analogue of the renormaliza- 
tion group equations, describing the cascade of gluon emission leading, to a 
power-law enhancement of soft gluons similar to the BFKL equationEZi The 
RLSM approach also describes screening effects in the parton distribution at 
small k±^ . The precise origin of this saturation has recently been elucidated by 
Kovchegovca 

The picture that emerges is the following: Gluons in the classical field 
generated by the valence quarks are fully Lorentz contracted by the Lorentz 
factor 7 associated with the colliding nuclei, but gluons spawned by a splitting 
of those primary gluons experience only a partial Lorentz contraction of order 
0:7, where x is the momentum fraction carried by the parent gluon. As the 
branching process evolves to softer and softer gluons, the spatial extent of this 
gluon cloud becomes more and more diffuse in the light-cone variable x^ . This 
result confirms the intuitive picture embodied in the original parton cascade 
model, and provides a quantitative formulation of it. 

4 Inmedium effects 

In-medium effects on parton-parton interactions are essential to the viabil- 
ity of the parton cascade model. The application of perturbative QCD to 
nucleon-nucleon collisions requires the introduction of ad hoc cut-offs describ- 
ing nonperturbative QCD effects, such as quark confinement and chiral sym- 
metry breaking. In-medium effects, which grow rapidly in size as function 
of A, produce perturbative cut-offs when the density of the medium becomes 
sufhciently high. For example, QCD is known to become perturbativeR at high 
temperature when the color-electric screening mass /id ^ Aqcd- 



"^Soine nonperturbative effects remain even at higti T, precisefy because static magnetic 
interactions are not screened by perturbative in-medium interactions. 



Although in-medium effects work in favor of the parton cascade model, 
they are not easily incorporated in practice. The problem is that in-medium 
effects are quite complicated and not easily treated correctly. The two main 
in-medium effects that are known to provide effective infrared cut-offs to per- 
turbatively divergent parton interactions are: 



• 



• 



color-electric screening, which suppresses soft 2-^2 scattering ampli- 
tudes; 

gluon radiation suppression, which reduces 2 — > 3 (and 2 ^ n) ampli- 
tudes with soft gluons in the final state. 



Dynamical screening, at lowest order, is described by the in-medium contri- 
butions to the one-loop gluon polarization function: 




At moderately high q±, the gluon population grows like n{k) « (^1^2)^/^, 
providing a screening scale ij-{A) that increases rapidly with the size of the 
nuclei participating in the coUisiontll 

Radiation suppression, also known as the Landau-Pomeranchuk-Migdal 
(LPM) effect, is a much more complicated mechanism. Its theoretical descrip- 
tion requires a good understanding of the multiple scattering problem in QCD. 
Considerable progress has recently been made in this area, especially through 
the work of Baier et alE3, ZakharovLJ and othersLJ The main difference be- 
tween QCD and the well-known case of QED is that a radiated gluon also 
rescatters in the medium at the same order in a^ as the radiating particle; this 
is not so for a photon radiated by a fast charge moving in an electromagnetic 
plasma. 

Revisiting the diagram shown in Fig. 0, the parton cascade model requires 
infrared cut-offs for both the central 2^2 scattering matrix element and each 
of the four branching cascades. This is where the in-medium effects help: at 
high density one expects that the modifications of the elementary scattering 
amplitude ensure an infrared safe behavior. To date, two attempts have been 
made to practically implement the action of these in-medium effects: 

1. In the self-screened parton cascade model (SSPCMo), the color-electric 
screening scale ^{pt) was calculated self-consistently for primary parton 
interactions only, and the further evolution of the parton plasma was 
described in the framework of the hydrodynamical model. 



2. In Klaus' PCM code VNItJ, the space-time picture of parton interactions 
is linked to the virtuality evolution of partons. A new interaction is 
permitted if its momentum transfer exceeds the virtuality scale of the 
participating partons at that time. Soft radiation is suppressed in the 
medium by this rule, because intermediate collisions continue to reset 
the parton virtuality to that of the latest collision. 

Both of these approaches are based on specific assumptions about the relation 
between the space-time and virtuality evolution of off-shell components of the 
parton distributions. One way to study this issue rigorously is the Wigner func- 
tion representation. First results c3 obtained by this method are interesting, 
but do not yet fully address the complications encountered in a QCD parton 
cascade. A more direct approach to the problem of the space-time picture of 
off-shell quantum fluctuation is based on a modification of the QCD evolution 
equations to include an infrared scaleoEZI In free space this infrared scale is 
determined by properties of the final state (hadronization scale) ; in a medium 
it is determined by screening effects. 

4-1 Self-screened parton cascade 

Here one considers the scattering of an initial state parton as completed after 
a time t{px) which depends on the momentum transfer in the reaction. The 
uncertainty relation suggests t{pt) ~ ^/pt- (We will drop the factor Ti in the 
following.) The scattered partons are then assumed to screen the scattering 
processes that involve a smaller momentum transfer: 

A^d(pt) = —o^sipl) / d'fc|V,n(fc)|. (4) 

The density of partons scattered at pr is, in turn, infiuenced by the screening 
because the differential cross section depends on /io: 

dpi {pi + fii{pT)y I ^ ' ^1 

If ^j-bIpt) becomes large enough at low px, so that da/dp^, remains perturba- 
tively small, the coupled set of equations can be integrated down to px = 0. 
Since the rapidity density of scattered partons grows as (^1^2)^/"^ Iny^s, this 
condition requires large A and high energy. The SSPCM concent-,has recently 
been investigated more formally by Makhlin and Surdutovicll53 within the 
framework of the closed-time-path formalism for real-time Green functions. 
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In this framework, the final state mass-shell condition for emitted partons reg- 
ulates the infrared divergences of single-particle correlation functions, such as 
the parton densities. Properly carried through, this concept leads to a self- 
consistent equation for the gluon screening mass (plasmon mass) similar to the 
one used in the SSPCM. 

Quantitatively, one finds that /ud approaches about 1 GeV at low px in 
Au -t- Au collisions at RHIC energy (100 GeV/u) and 1.5 GeV at LHC energy 
(2.75 TeV/u). The differential minijet cross section as function of pr peaks 
at about the same value, clearly showing the improved infrared behavior of 
the self-screened parton cascade. The total deposited energy within one unit 
of rapidity and after a characteristic formation time of 0.25 fm/c is eo ~ 60 
GeV/fm3 (RHIC) and eo ~ 430 GeV/fm^ (LHC). The conditions established 
by the SSPCM can, therefore, be taken as initial conditions for the thermal 
and chemical evolution of a quasi-equilibrated parton plasma. 

The kinetic equations for the evolution of such a plasma were-jderived by 
Biro et al. 123 and by Xiong and Shuryaktj. Extensive calculationsc2l, including 
longitudinal and transverse expansion, have shown that the plasma cools down 
to the critical temperature of QCD (T^ ~ 150 MeV) after 5 fm/c (RHIC) and 
10 fm/c (LHC). The emission of electmmagnetic probes by such an evolving 
QCD plasma has also been calculatedc3 

4-2 Monte-Carlo space-time cascade 

The statistical implementation of parton cascades in the Monte-Carlo simu- 
lation code VNIE3 achieves an improved infrared behavior through heuristic 
rules that suppress certain interactions on the basis of kinematic considera- 
tions. The first rule asserts that independent scattering events involving the 
same parton require a sufficient time separation so that the time between scat- 
terings is larger than the duration of the individual events. With the duration 
of an interaction again defined as t{pt) ^ p^ , where pT is the momentum 
exchange, this requires that the time between interactions At > t{pt). An- 
other way of ensuring this condition is to endow a parton after a scattering 
by pt with an initial virtuality Qo = Pt, which then gradually decreases with 
time as Q{t) = Qot{pt)/t. A subsequent scattering with pip requires that 
pip > Q{t) at the moment of the interaction. A second similar rule suppresses 
soft parton splittings in the presence of multiple scattering. Again, this rule 
can be formulated in terms of a parton virtuality that decreases with time 
between scatterings and is reset by each new interaction (see Fig. ^. 

Although VNI still contains "arbitrary" infrared cut-off parameters (deter- 
mined by a comparison with nucleon-nucleon interactions), these are needed to 
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/o/ \ ^ Figure 4: Illustration of in-medium sup- 

^:l' J q ^T g^^'^ ;. pression effects incorporated in the VNI 

3 ^^^ code. The virtuality of a scattered parton 

Pj- 9 evolves with time, Q(t). Sequential scatter- 

ed ings or branchings are suppressed, if Q{t) 

is too large. 



limit soft scatterings in the initial set of parton interactions, and they become 
important toward the end of the cascade evolution when the parton plasma 
becomes more and more dilute. The cut-off parameters effectively determine 
the end of the cascade evolution and, by suppressing soft interactions, they 
trigger the hadronization of parton clusters. 

The fact that this improved parton cascade is much less dependent on arbi- 
trary cut-off parameters makes it possible to explore applications of this model 
at lower energies where nucleon=iiucleon collisions provide little information. 
Recently, Geiger and SrivastavaciJ have studied the predictions that the VNI 
code makes for nuclear collisions in the energy regime of the CERN-SPS. While 
most of the particle yield at rapidities |y| > 2 is produced by fragmentation of 
the unscattered beam remnants, the model predicts a significant contribution 
to particle production a± central rapidity from partons that have undergone 
perturbative scatteringu3 This contribution is rising rapidly with nuclear mass 
A, roughly as (^1^2)^^'^. This perturbative contribution to the energy de- 
position at \y\ < 1 coincides with a rapid increase of the energy density in 
scattered partons at r < 1 fm/c, which rises from about 2 GeV/fm^ in S-l-S to 
5 GeV/fm'^ in Pb -I- Pb. This rise may be correlated with the much enhanced 
suppression of charmaaium production in Pb + Pb collisions as observed by 
the NA50 experimente^l 

5 Open problems 

In spite of the considerable conceptual advances in the parton cascade model 
compared to its original formulation, certain aspects have still not been com- 
pletely clarified. For issues concerning the initial state I can only refer to the 
recent review article of McLerran and VenugopalanCj One problem associated 
with the final state concerns the source of the produced entropy: The classical 
gluon field radiated in the encounter of two specified color chaxee distributions 
retains its full coherence during the course of the interactioncJ The Gaussian 
average over inital conditions formally introduces entropy at late times, but a 
significant part of this entropy is already present in the initial state. The prob- 
lem of thermalization remains a theoretical challenge even in this framework. 
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It is likely that its solution must be sought in the chaotic dynamical aspects of 
classical Yang-Mills fields that have beea-Qbserved in numerical simulations of 
the classical non-abelian gauge theorycZio The apparent decoherence of clas- 
sical gauge fields is visible in a dramatic fashion when the collision of t. 
non-abelian wave packets is studied in numerical simulations on the latticel 

A consistent implementation of the RLSM into the parton cascade model 
requires a transport description of partonsJxiLluding mean color fields. Such a 
framework has been known for a long timetjO, but practical implementations 
of the QCD Bpltzmann equation for partons and fields have been attempted 
only recentlytj The idea behind this approach is to separate short-distance and 
long-distance dynamics by means of a lattice cut-off: excitations with momenta 
k < n/a, where a is the lattice spacing, are represented as classical fields on 
the lattice; those with higher momenta are represented as colored particles. 
The lattice cut-off must then be chosen such that g^ < ir/a < fi. "Hard" 
collisions, i.e. those with a momentum transfer q > 7r/a, are now infrared safe 
due to the lattice cut-off, soft collisions are represented as interactions among 
particles via the lattice fields. These incorporate the crucial screening effects 
at the scale g/i. 

Another conceptual issue that is understood, in principle, but whose de- 
tailed investigation is an important outstanding problem, is the question what 
constitutes a consistent set of semiclassical transport equations for a theory 
with perturbatively massless modes such as QCD. This issue is conceptually 
resolved in the case of a massless scalar quantum field theory, where the in- 
troduction of a self-consistent medium-dependent mass term is sufficient. The 
issue is trickier in the case of QCD because of two problems: the need to retain 
gauge invariance and the long radiative tails of the spectral functions for the 
colored quasi-particles. Gauge invariance essentially requires that every mod- 
ification of a n-body correlator is mirrored by an analougous modification of 
the (n -I- l)-body correlator satisfing Ward idenities. The treatment of radia- 
tive tails of the spectral function, on the other hand, requires renormalization 
group techniques. 

At the time of his death, Klaus was working on this problem, applying 
the "exact" renormalization group technique Ej to QCD transport theory. On 
the evening before he boarded iJie ill-fated Swissair flight, he sent an almost 
finished draft of a manuscript EJ to Wetterich and me with the request: "I 
expect to see your comments when I return." I believe that, again, Klaus was 
onto something very important, but it will remain for others to finish what 
he started. Let me just end this part of the discussion with the remark that 
the usual mass-shell condition in transport theory is intimately connected with 
the renormalization-group equations, and that a better understanding of this 
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connection is crucial to further progress in the quantum transport theory of 
QCD. 

The next logical step in any parton cascade model that aims at practical 
predictions for experimentally observables is hadronization. InjxtUaboration 
with John Ellis, Klaus worked intensely on this difficult problemLj Their work 
is based on an effective theory of scale invariance breaking in QCD and has 
had some success in describing hadronization of quark jets produced in e+e~ 
annihilations. Further progress here will have to rely on an improved under- 
standing of the mechanism of quark confinement in QCD. 

Finally, let me point out that Steffen Bass had been in the midst of a coUab- 
pjjation with Klaus to add a state-of-the-art hadronic cascade model (UrQMD 
Ea) to the parton cascade code VNI, in order to describe final-state interactions 
of the hadrons created by hadronization of the quark-gluon plasma. The first 
results of this project are intriguing: they indicate that a small fraction of the 
created hadrons emerge without rescattering, whereas the bulk ofj-the hadron 
yield is reprocessed through the melting pot of hadronic reactionsEH 

6 Summary 

The parton cascade model was developed by Klaus to provide a QCD-based 
description of the approach to a locally thermalized state in collisions of heavy 
ions in the RHIC energy regime and beyond. In its original formulation the 
PCM predictions were critically dependent on several cut-off parameters that 
had to be determined from pp collision data. Recent advances incorporating 
in-medium effects into the parton interactions have reduced this dependence 
significantly, possibly allowing the application of the PCM over a wider en- 
ergy range. Results obtained for nuclear collisions at CERN-SPS energies are 
intriguing. 

The in-medium effects that modify parton-parton interactions not only 
reduce the parameter dependence of the model, they also provide valuable in- 
sight into the dynamics of a dense parton plasma. It is clear that we are just 
at the beginning here. The transport properties of off-shell quanta need to be 
understood much better, not only in cases where the off-shell propagator is 
dominated by a well-defined resonance, but especially in the case where the 
particles never get close to their mass shell as it applies to QCD. Another open 
question concerns the need for mean color fields. Such fields are not included 
in present versions of the PCM, but the random light-cone source model sug- 
gest that mean fields may be essential ingredients of a complete description 
of soft processes in nuclear collisions. One needs to take an average over a 
Gaussian ensemble of mean fields, where the width of the field distribution 
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is more important than the expectation value which remains zero. It would 
be interesting to explore possihk. pjnnections of the RLSM to the traditional 
chromo-hydrodynamical modelcJ'cEl 

Ultimately, the question is whether the parton cascade model can be re- 
placed by a controlled approximation scheme where, in principle, successive 
orders of ever more sophisticated corrections can be calculated. We are still 
some steps away from a consistent formulation of transport phenomena off 
equilibrium in QCD. It is even unclear whether we even know what the small 
parameters in such an approximation scheme are. It is clear that a high den- 
sity of excitations of the QCD vacuum is an essential condition, but there are 
many subtleties if one wants to go beyond this statement. However, steady 
progress in this field is being made, and there is reason to hope that a consistent 
formulation of transport phenomena in QCD can ultimately be achieved. 

Finally, the treatment of the late phase of a rclativistic heavy ion collision, 
when the dense matter breaks up into individual hadrons, has recently seen 
some exciting improvements. One still needs to rely on a phenomenological 
hadronization model for the transition to a hadronic cascade, but the latter 
can be studied with he technology that has been developed for the description 
of nuclear collisions at lower energies. In fact, the description of the hadronic 
final-state interactions should be much more reliable than, e.g. that of a nuclear 
collision at AGS energies, because the hadronic system starts out close to 
thermal equilibrium. Klaus was very excited about this approach, which will 
allow for a better comparison of the predictions of the parton cascade model 
with the experimental data that will soon become available from RHIC. 
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